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Abstract. We discuss rare meson decays K + — > ir~£ + £' + and D + — » K~i Jr l' Jr 
(£, t! = e, /i) in a supcrsymmetric extension of the Standard Model with i?-parity 
violation. Estimates of the branching ratios for these decays are presented. 

In the standard model (SM), lepton L and baryon B number conservation 
take place due to the accidental J7(1)l X U(1)b symmetry existing at the level 
of renormalizable operators. But, for many extensions of the SM, this is not the 
case. The well known mechanism of lepton number (LN) violation is based on 
the mixing of massive Majorana neutrinos predicted by various Grand unified 
theories (GUTs) [1]. The Majorana mass term violates LN by AL = ±2 [2] 
and can lead to a large number of LN violating processes. Among them, the 
most sensitive to the LN violation are processes such as neutrinoless double 
beta decay (A, Z) — > (A, Z + 2) + e~ + er (for a recent review, see [3]), rare 
meson decays (see, e.g., [4,5]) 

M+ -> m'-£ + £+ , (1) 

and like-sign dilepton production in high-energy hadron-hadron and lepton- 
hadron collisions (see, e.g., the papers and references therein: pp — > i^l'^X [4, 
6], e + p — > i? e £ + £' + X [7]), which have been extensively studied in the literature. 

There exists now convincing evidence for oscillations of solar, atmospheric, 
reactor, and accelerator neutrinos [8]. The oscillations, i.e., periodic neutrino 
flavor changes, imply that neutrinos have nonzero masses and they mix among 
each other: neutrinos vt of specific flavors I = e, fx, r are the coherent superpo- 
sition of the neutrino mass eigen-states vn of masses toat, 

vt = Uinvn- (2) 

Here the coefficients Uin are elements of the unitary leptonic mixing matrix - 
the PMNS matrix [9]. 

Neutrino flavor changes imply lepton family number Lg nonconservation ad- 
missible for neutrinos of both types, Dirac and Majorana, but for the Dirac neu- 
trinos, in contrast to the Majorana ones, the total lepton number L = Lg is 
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conserved. The nature of the neutrino mass is one of the main unsolved prob- 
lems in particle physics. However, oscillation experiments can not distinguish 
between the two types of neutrinos. 

In Refs. [4,5] we investigated the rare decays (1) of the pseudoscalar mesons 
M = K, D, D Sl B, mediated by light (m^r <C me, my ) and heavy (m^v 3> Tom) 
Majorana neutrinos. It was shown that the present direct experimental bounds 
on the branching ratios (BRs) are too weak to set reasonable limits on the 
effective Majorana masses. Taking into account the limits on lepton mixing 
and neutrino masses obtained from the precision electroweak measurements, 
neutrino oscillations, cosmological data and searches of the neutrinoless double 
beta decay, we have derived the indirect upper bounds on the BRs that are 
greatly more stringent than the direct ones. 

In this report, we investigate another mechanism of the AL = 2 rare decays 
(1) based on i?-parity violating supersymmetry (SUSY) (for a review see [10]). 
We recall that i?-parity is defined as R = (— 1) 3 ( B_L )+ 2S , where S, L, and B are 
the spin, the lepton and baryon numbers, respectively. In the minimal super- 
symmetric extension of the SM (MSSM) , i?-parity conservation is imposed to 
prevent the L and B violation; it also leads to the production of superpartners 
in pairs and ensures the stability of the lightest superparticle. However, neither 
gauge invariance nor supersymmetry require i?-parity conservation. There are 
many generalizations of the MSSM with explicitly or spontaneously broken R- 
symmetry [10]. We consider a SUSY theory with the minimal particle content 
of the MSSM and explicit imparity violation (#MSSM). 

The most general form for the i?-parity and lepton number violating part of 
the superpotential is given by [10, 11] 

W n = e afj (±\ ijk L?L?E k + \' ijk LfQ^D k + aLfH^j . (3) 

Here k — 1,2,3 are generation indices, L and Q are SU(2) doublets of 
left-handed lepton and quark superfields (a,f3 = 1,2 are isospinor indices), 
E and D are singlets of right-handed superfields of leptons and down quarks, 
respectively; H u is a doublet Higgs superfield (with hypercharge Y = 1); \ij k = 
—Xjik, X'^f, and ej are constants. 

In the superpotential (3) the trilinear (oc A, A') and bilinear (oc e) terms 
are present. In this work, we assume that the bilinear terms are absent at 
tree level (e = 0). They will be generated by quantum corrections [10], but it 
is expected that the phenomenology will still be dominated by the tree-level 
trilinear terms. 

At first we consider the rare decay K + (P) — > tt~(P ) + £ + (p) + £ + {p ) in 
the $MSSM (a rough estimate of the width of the decay B(K+ -> 7r"/i+^+) 
in the same theory was obtained in [12]). The leading order amplitude of the 
process is described by three types of diagrams shown in Fig. 1. 

The hadronic parts of the decay amplitude are calculated with the use of a 



2 



f\ If 




Figure 1: Feynman diagrams for the decay K + —* n~ +l + +1 + mediated by Majorana neu- 
trinos v, neutralinos x°, gluinos g with / being the scalar superpartners of the corresponding 
fcrmions / = £,u,d (leptons and quarks). Bold vertices correspond to Bethe-Salpeter am- 
plitudes for mesons as bound states of a quark and an antiquark; (t), (b) and (3) stand 
for tree, box and simply the third kind of diagrams, respectively. There are also crossed 
diagrams with interchanged lepton lines. 



model for the Bethe-Salpeter amplitudes for mesons [13], 

Xp{q)^l b {l-5 M rWp{q), (4) 

where Sm = (mi + m 2 )/m M , tum is the mass of the meson composed of a 
quark q\ and an antiquark q-i with the current masses mi and 777,2, P = Pi + Vi 
is the total 4- momentum of the meson, q = (jpi —p2)/2 is the quark-antiquark 
relative 4-momentum; fp(q) is the model-dependent scalar function. 

For all mesons in question, vtlm <C msusY, where msuSY ^ 100 GeV is 
the common mass scale of superpartners, and for heavy Majorana neutrinos, 
777 jv ^ 777 jvf (the contribution of light neutrinos is strongly suppressed by phe- 
nomenology [4,5]), we can neglect momentum dependence in the propagators 
(see Fig. 1) and use the effective low-energy current-current interaction. In 
this approximation the decay amplitude does not depend on the specific form 
of the functions (fip(q) (see Eq. (4)) and is expressed through the known decay 
constants of the mesons, /m, as 



fit = 4 V^c 5 M (2tt) 4 j d 4 q <p P (q) 



where N c — 3 is the number of colors. 

For the total width of the decay we obtain 
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Here $«< is the reduced phase space integral (2 = (P — P') 2 /m 2 K ): 

<&«, - ^ ~ (l - l(h+ - z){h- - z)(l + - z){l_ - z)} 1/2 , 

and the various parameters are defined as follows: 

h± = (1 ± m^/m K ) 2 , Z± = [(me ± m^/m^] 2 ; 
€lsW = -T^)N S2 + tanMW) - Q(il>))N S i, 
emiip) = Q(ip) tan OwNsi, 

where Q{ip) and T 3 (i/>) are the electric charge and the third component of 
the weak isospin for the field ip, respectively, and Ns a is the 4x4 neutralino 
mixing matrix. For the numerical estimates of the branching ratios, B«/ = 
T (M + — > M'~£ + £' + ) /r to tai, we have used the known values for the couplings, 
decay constants, meson, lepton and current quark masses [5,8], and a typical set 
of the matrix elements N$i,Ns2 from Rcf. [14]. In addition, we have taken all 
the masses of superpartners to be equal with a common value msusY ■ Taking 
into account the present bounds on the effective inverse Majorana masses [5], 
we find that the main contribution to the decay width comes from the exchange 
by ncutralinos and gluinos (see Fig. 1). The results of the calculations with 
the use of Eq. (5) for the decays K + — ► ir~£ + £ + , and an analogous formula 
for the decays D + — > K~£ + £' + , are shown in the fourth column of Table 1 
(here TO200 = ws(/sy/(200GeV)). In the second and third columns of this 
table, the present direct experimental upper bounds on the BRs [8] and the 
indirect bounds for the Majorana mechanism of the rare decays [5] are shown, 
respectively. Our result for the B(K + — > 7r~/i + /i + ) is in agreement with a 
rough estimate of Ref. [12]. 

To calculate the upper bounds on the BRs in the $MSSM, we take m 2 oo = 1 
and|A^ fe A^., fc ,|<10- 3 [15]. It yields 

B (K+ - ■*-£+£'+) < 1(T 23 , B (D+ -> K-£+£'+) < KT 24 . 

These estimates are much smaller than the corresponding direct experimental 
bounds but are close (except for the ee decay mode) to the indirect bounds 
based on the Majorana mechanism of the decays (see Table 1). 
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Table 1: The branching ratios B , for the rare meson decays M+ — > M' . 



Rare decay 


Exp. upper 
bound on B«< 


Ind. bound 
on B«> M SM) 


B«' x m^ 
(/?MSSM) 


K+ -^n-e+e+ 


6.4 x 10" 10 


5.9 x 10~ 32 


5.7 x 10- iY |A' ni A' 112 | 2 


K + — > ir~ fi + /i + 


3.0 x l(T y 


1.1 x 10~ 24 


2.0 x 10- iY |A 211 A 212 2 


K + — > 7r~e + /i + 


5.0 x 10" 10 


5.1 x 10~ 24 


1.9xl0^|Ai n A 212 +A 211 Ai 12 | 2 


L>+ -» if-e+e+ 


1.2 x 10~ 4 


1.5 x lO" 31 


1.0 x 10- iS |Ai 22 Ai n + 0.45A' 121 A' 112 | 2 


D+ -> K- 


1.3 x 10" b 


8.9 x 10~ 24 


9.6 x 10- iy |A 222 A 211 +0.45A 221 A 212 | 2 


D+ -» ^-6+^+ 


1.3 x 10~ 4 


2.1 x 10~ 23 


4.9xl0-H(Ai 22 A 211 +A' 222 A' l2ll ) 
+0.45(A 121 A 212 + A 221 A 112 ) 
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